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The generalized Rabi model is exactly solved by employing the unitary transformation method in 
the occupation number representation. The analytical expressions for the complete energy spectrum 
consisting of two sub-energy spectra are presented in the whole range of all the physical parameters. 
Each energy level possesses twofold degeneracy and is determined by the parameter in the unitary 
transformation, which obeys a highly nonlinear equation. The corresponding eigenfunction is a 
convergent infinite series in terms of these physical parameters. Due to the level crossings between 
the neighboring eigenstates at certain parameter values, such the degeneracies could lead to novel 
physical phenomena in these quantum systems with the light-matter interaction. 
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The Rabi model describes a two-level system coupled 
with a single bosonic mode [1]. Such a simplest interact- 
ing quantum model has had wide applications in many 
fields of physics, e.g. atomic physics [2], quantum optics 
[3], trapped ions [4,5], quantum dots [6], superconducting 
qubits [7,8,9], cold atoms [10], and etc.. It is also expect- 
ed to be the theoretical basis for quantum information 
and quantum technology [11-14]. 

A generalized Rabi model usually has the Hamiltonian 
(15) 


H = wata + gla +a)o, + Aoz + 04, (1) 


where o, and g, are the Pauli matrices for the two-level 
system with level splitting 2\, at and a are the creation 
and annihilation operators for the single bosonic mod- 
e with frequency w, respectively, the light-matter inter- 
action is controlled by the coupling parameter g, and 
the last term eo, is the driving term which breaks the 
Z2 symmetry and leads to tunnelling between the two 
levels. We note that the competition between g and w 
produces the different experimental regimes. When g/w 
is small, by applying the rotating-wave approximation, 
the Rabi model (1) with e = 0 leads to the so-called 
Jaynes-Cummings model [16], which is relevant to most 
experimental regimes. Because the Jaynes-Cummings 
model is integrable, it is easy to derive its analytical so- 
lution. With increasing g/w, the ultrastrong coupling 
regime (g/w >~ 0.1) [12] or the deep strong coupling 
regime (g/w >~ 1.0) [9] is reached, where the Jaynes- 
Cummings model is invalid and cannot be used to inves- 
tigate the interaction between light and matter. Recent- 
ly these regimes have rapidly growing interesting due to 
their fundamental characteristics and the potential ap- 
plications in quantum devices [11-14]. 

Although the Hamiltonian (1) has a simple form, it 
has not been possible to obtain its analytical solution, 
which is considerably important for exploring accurately 
the light-matter interaction, from weak to extreme strong 
interaction. In Ref. [15], Braak presented an analytical 
solution of the Rabi model by using the representation 


of bosonic operators in the Bargmann space of analytical 
functions. The energy spectrum consists of two parts, the 
regular and the exceptional spectrum. However, such a 
spectrum structure is considerably strange. It has been 
proved that the Braak’s analytical solution of the Rabi 
model is completely wrong due to the derivation error in 
solving the time-independent Schrodinger equation [17]. 

In this work, we exactly diagonalize the Hamiltonian 
(1) in the whole range of the physical parameters by us- 
ing the unitary transformation technique in the occupa- 
tion number representation. Such a direct and powerful 
approach has been used to solve successfully the com- 
plex two-dimensional electron gas in the presence of both 
Rashba and Dresselhaus spin-orbit interactions under a 
perpendicular magnetic field [18]. The eigenvalues are 
fixed by the parameters in the unitary transformations, 
which satisfy highly nonlinear equations. The eigenfunc- 
tions are convergent infinite series in terms of the physical 
parameters. The exact Landau levels have been accepted 
popularly [19](also see the citations of Ref. [18]). When 
the Rashba or Dresselhaus spin-orbit coupling vanishes, 
this two-dimensional electron system has the same boson- 
ic Hamiltonian with the Jaynes-Cummings model except 
a constant [16,18]. 

The two-component eigenstate of the Hamiltonian (1) 
for the nth energy level with quantum number s usually 
has the general form 


Mu = (L+ AR.) Do (ome? + 18m 77), (2) 


m=0 


where the 2 x 2 matrix is a unitary one, s = +1, Ans 
is a real parameter to be determined below by requiring 
ars and 875 to be nonzero, dm is the eigenstate of the 
mth energy level in the occupation number representa- 


tion, i.e. atom = VM + lọm+1, Adm = VMbm_—1 and 


< bm'|ém >= mm. When m —> +00, at = p25 = 0. 
We solve the eigen-equation H|n,s >= Ens|n,s > and 
obtain the double degenerate energy spectrum consisting 
of two sub-energy spectra, which is very different from 
that presented in Ref. [15]. 

The sub-energy spectrum I. In order to obtain the an- 
alytical solution of the Hamiltonian (1) in the whole pa- 
rameter space, we first choose 


—A?2 
[u(n +1) + rate Tew €— Ensag 
reaetovn FIER = 0, 


zl Ens) 873 
+ {cae rga + lan, = 0 


in the eigen-equation. Then the eigenfunction associated 
with the eigenvalue Ens is solely fixed by requiring 


PAS + (1 — A2,)]82° — 2gAnsVn + lars, =0, (4) 


or 


[2AAns T e(1 = A? ,)Jar and F 2gAns vn + 13, = = 0. (5 ) 


We solve the homogenous linear equations (3) about 
ans, and 6p? by vanishing of the coefficient determi- 
nant. Then the eigenvalue for the nth eigenstate with s 
has the analytical expression 


Ens = (n+ $)w + 8=Ens, 


Note that the quasiparticle energy Ens must be larger 
than zero. Combining Eqs. (3)-(5), the parameter A,,, 
is determined by the constraint 
e(1 + A?,) — 2Ans(Ens — wn) = 0, (7) 
or 
e(1 + A?,) +2Ans[Ens — w(n + 1)] = 0. (8) 


Surprisingly, Eq. (7) with s = —1(1) coincides with Eq. 
(8) with s = 1(—1). So we have 


e(1 T Ahs) A Ans (20Ens = 


where ø = +1. Obviously, when e = 0, then Ans = 0 
from Eq. (9). Therefore, the eigenvalue (6) has a simple 
form in the absence of the driving term. Fig. 1 plots the 
low-lying energy levels as a function of g at different A 
under € = 0. We can see that there are level crossings 
between the neighboring eigenstates. With increasing A, 
the energy levels with s = 1(—1) become higher (lower), 
and these crossing points move toward the origin. 
When e€ # 0, Ans in Eq. (9) with o = 1 has an w- 
dependent solution. The corresponding eigenvalues Ens 


(a) A/o=0.3 


FIG. 1: (Color online) The low-lying energy levels of the sub- 
energy spectrum I in unit of w as a function of the coupling 
parameter g at different A under e = 0. The solid lines denote 
n = 0,1,--- ,5 and s = 1 while the dash lines mean n = 
1,2,--- ,5 and s = —1. 


(n = 0,1,2,--- ,co,s = +1) form the sub-energy spec- 
trum I. In Fig. 2, we show the low-lying energy levels 
of the sub-energy spectrum I as a function of g at dif- 
ferent À under e = 0.4 and the corresponding parameter 
Ans. From Fig. 1 and Fig. 2, we find that if A or e 
increases, the energy levels with s = 1(—1) also move 
up (down), and the crossing points also shift toward the 
left side. When g = 0, the sub-energy spectrum I (6) 
recovers the exact eigenvalues for the interactionless case 
Ens = wn+sVd2 +e with Ans = (A — VX + €)/e. 
Obviously, if € > 0, Ans > 0. We note that another 
sub-energy spectrum corresponding to the w-dependent 
solution Ans in Eq. (9) with o = —1 is nothing but the 
sub-energy spectrum II (see below). 


For the eigenstate associated with the eigenvalue (6), 
we have 


ns __ (1 + A?) (Ens — nw) + (1 — A? JA — 2Ans€ ns 


An+1 = vn+1(1 = A2.)g n? 

(10) 
where (7'* is an arbitrary constant and can be set to 1, 
and the coefficients a?* and 675 are uniquely determined 
by the recursion relations 


ive) = (E vero E) 
ii a 


for i = 0,1,2,- 


VIFM” ( i ) = syne ( FY 
my E Ag 


n, and 


i—1 
(12) 
for i =n+1,n+2,--- ,+00. Here we have defined 
Ms = rraz (1 E Ae oy = 2AnsOzl, 
ns 4 i ns = Ans 
N; = (wi Ens )I T Conte A4 T 1+42, €)Ox (13) 


1-42, 2A, 
+(qa8A -= IFAS, €)Oz, 


where I is the 2 x 2 unit matrix. 
The sub-energy spectrum II. Now we use another choice 


glo glo 


FIG. 2: (Color online) The low-lying energy levels of the sub- 
energy spectrum I in unit of w as a function of the coupling 
parameter g at different \ under e€ = 0.4, shown in (a) and (c). 
The solid lines denote n = 0,1,- -- ,5 and s = 1 while the dash 
lines mean n = 1,2,--- ,5 and s = —1. The corresponding 
Ans are displayed in (b) and (d), respectively. 


-A2 
(wn + IFR A E e — Ens)an? 
ns I-A? ns 
; +izartgvn + 1 mi = 0, (14) 
-A 
[w(n + 1) tat’ H ewes €— Ens] Bni 


+a gvn + lags = 0 
in the eigen-equation. The corresponding eigenfunction 
is uniquely determined by letting 
[2AAns +e(1— A? )lar® +2gAn./n +1824; =0, (15) 
or 
[2>Ans +€(1— A2,)]bn4 , —2gAnsVn + 1a% =0, (16) 


Solving Eqs. (14), we have the eigenvalue for the nth 
eigenstate with s 


Ens = (n+ $)w+ Ons, 


w 1 1-A?, ia 1—42; 
Ons = y (9 + IRRA — Pare)? + (n + 1)9 (FA)? 
(17) 
Here Ans satisfies the nonlinear equation 
e(1 + A?) + 2Ans(Ens — wn) = 0, (18) 
or 
e(1+A?2,) —2Ans[Ena —w(n +1) =0, (19) 


which is derived from Eqs. (14)-(16). Similar to Eqs. (7) 
and (8), Eq. (18) with s = —1(1) is also consistent with 
Eq. (19) with s = 1(—1). Therefore, we obtain 


e(1 + A?) + Ans(2TOns +w) = 0, (20) 


(b) Aa=0.5 
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FIG. 3: (Color online) The low-lying energy levels of the sub- 
energy spectrum II in unit of w as a function of the coupling 
parameter g at different A under e = 0. The solid lines denote 
n = 0,1,--- ,5 and s = 1 while the dash lines mean n = 
0,1,---,5 and s = —1. 


where 7 = +1. If e = 0, then Ans = 0 in Eq. (20). So 
the eigenvalue (17) also has an explicit expression in this 
case. We depict the low-lying energy levels as a function 
of g at different \ under e = 0 in Fig. 3. With increasing 
À, the energy levels with s = —1(1) become higher (lower) 
while the crossing points move away the origin, opposite 
to those in Fig. 1. 

If €e £0, Ans in Eq. (20) with 7 = 1 also has an w- 
dependent solution. The corresponding eigenvalues con- 
stitute the sub-energy spectrum II. Fig. 4 exhibits the 
low-lying energy levels of the sub-energy spectrum II as 
a function of g at different \ under e€ = 0.4 and the cor- 
responding parameter Ans. When A is increasing, the 
energy levels with s = 1(—1) also move up (down), and 
the crossing points also shift toward the origin. This re- 
sembles the level shifts in the sub-energy spectrum I (see 
Fig. 2). From Figs. 3 and 4, we find that when à = 0.3w, 
the energy levels with s = —1(1) move up (down) with in- 
creasing €. Oppositely, the energy levels with s = 1(—1) 
move up (down) at A = 0.5w. Obviously, the level shift- 
s in the sub-energy spectrum II are more complex than 
those in the sub-energy spectrum I. We note that the 
sub-energy spectrum associated with Ans in Eq. (20) 
with 7 = —1 is consistent with the sub-energy spectrum 
I. Therefore, both the sub-energy spectrum I and II are 
double degenerate. So the Braak’s mathematical solution 
also loses a half of the eigenvalues. 

For the nth eigenstate with s in the sub-energy spec- 
trum II, we have 


ns (1 Eg A? (Ens a 


nw) — (1— A2,)A+ 2Anse 


ns 


= Q. ; 

= vn + 1(1— Aks) " 
(21) 
where a} is an arbitrary constant and is set to 1. The 


coefficients a? and 8?° obey the same recursion relations 
(11) and (12) in the sub-energy spectrum I. 

In summary, we have exactly solved the generalized 
Rabi model (1) in the framework of quantum mechanics. 
The complete spectrum is comprised of two double-fold 
degenerate sub-energy spectrum I and II. Very different 
from the energy spectrum in Ref. [15], the level crossings 
occur even if € is not a multiple of w/2. Such the physi- 
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FIG. 4: (Color online) The low-lying energy levels of the sub- 
energy spectrum II in unit of w as a function of the coupling 
parameter g at different À under e€ = 0.4, shown in (a) and (c). 
The solid lines denote n = 0,1,--- ,5 and s = 1 while the dash 
lines mean n = 0,1,--- ,5 and s = —1. The corresponding 
Ans are displayed in (b) and (d), respectively. 


cal solution can help us to deeply understand the light- 
matter interaction, especially in strong coupling regime. 
It is expected that the level crossings at certain physical 
parameters could produce novel physical phenomena in 
the interaction systems, just like in two-dimensional elec- 
tron gas with spin-orbit interaction under a perpendicu- 
lar magnetic field [20-22]. Finally, we would like to men- 
tion that the unitary transformation method presented 
here can be applied naturally to the general Rabi mod- 
el, which is a further extension of the generalized Rabi 
model by adding H’ = (at + a)(gyoy + gzoz) + €yoy to 
the Hamiltonian (1), where gy, gz, and ey are the physical 
parameters. 
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